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Abstract 

A novel technique is used to compute the bulk viscosity of high temperature holographic 
gauge theory plasmas with softly broken conformal symmetry. Working in a black hole 
background which corresponds to a non-trivial solution to the Navier-Stokes equation, and 
using a Ward identity for the trace of the stress-energy tensor, it is possible to obtain an 
analytic expression for the bulk viscosity. This can be used to verify the high temperature 
limit of a conjectured bound on the bulk viscosity for these theories. The bound is saturated 
when the conformal symmetry-breaking operator becomes marginal. 
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1 Introduction and summary 



The gauge-gravity duality [H [2j [3] has provided an invaluable tool for studying the strongly 
coupled limit of large N gauge theories. Perhaps surprisingly, there are features of these 
gauge theories which may be relevant to the quark-gluon-plasma (QGP) created at the Rel- 
ativistic Heavy Ion Collider (RHIC) at the Brookhaven National Laboratory. In particular, 
it seems that the shear viscosity to entropy density ratio of the QGP formed at RHIC is of 
the order of [I], 

2 =7- d) 
with rj the shear viscosity and s the entropy density. The result in ([1]) was first computed 
in [5] and was later shown to hold in a large class of gauge theories with a holographic 
dual [HI LTl IH] - Similar relations exist for the electrical conductivity [H] and for higher order 
transport coefficients [TO]. In a conformal gauge-theory plasma the bulk viscosity, £, vanishes, 
but in a non-conformal theory the bulk viscosity is, in general, non-zero. In [11] it was 
conjectured that the bulk viscosity always satisfies 

with d the number of space-time dimensions and c s the speed of sound, for large N, strongly 
coupled gauge theories with a holographic dual. The bound (j2J) was shown to be satisfied in 
[T2| [TT| [T3| [HI [T5J [16] for various configurations. In [17] an explicit setup where the bound is 
slightly violated was constructed. It is not clear if the bulk matter content of this particular 
example follows from a truncation of string theory. 

In this note, I show that in gauge theories with holographic duals involving gravity and a 
scalar which softly breaks the conformal symmetry, the bulk viscosity £ and speed of sound 
c s are given by 

cl = (d - 2A) tan ( ^) D(A, d) ( ±) ^ (l + o(±)) (3) 



where 



d-1 s y ' \ d ) y ' ' \T ) \ \T 

d \ d ~ 1 /A\ 2< - d_A ' / /A 

■£r) C = 2vr(rf-A) J D(A,rf)f-j I 1 + ° U ) 



16^-A)^-a-D / r(f) 

V ^ d >- 2(47T)^-A)( rf -l)2 l r (f + |)' ' 



1 



A is the scaling dimension of an operator dual to a bulk scalar, and A d ~ A is the source term for 
this operator which is what breaks conformal invariance. It is assumed that (d—2)/2 < A < d 
so the scalar operator is relevant and above the unitarity bound. An expression similar to 
([3]) was obtained in [181 EE] for d = 4 and 2 < A < 4 and expression (4J) restricted to d = 4 
and 2 < A < 4 was conjectured in [20], based on a numerical analysis] 1 ] Since (T4J) is valid in 
the high temperature limit, in theories with multiple scalars one should add a contribution 
of the form (jlj) for each scalar. The expressions for the bulk viscosity and speed of sound 
agree with the numerical results in the literature [22] ^1 

Using ([3]) and (jlj) and computing the leading contribution to the shear viscosity r], 



d 



d-l 



one finds that 



C/t] d-A /ttA 

1 rr" = -27TCOt — — 

r -~ d-2A \d 



lim 



> 2. 



(6) 



Curiously, the bound ([2]) is saturated in the limit where the conformal dimension of the 
scalar operator becomes marginal, A = d. 



2 Computation of the bulk viscosity 

In the hydrodynamic approximation the dynamics of a fluid are specified by a velocity field 
u^ix) (canonically normalized so that u^u^ = —1), an energy density e(x) and chemical 
potentials In what follows the chemical potentials will be set to zero. The energy 

momentum tensor of viscous hydrodynamics takes the form 

= eu"u v + PP» V - 7]a" u - (P^d a u a + 0(d 2 ) (7) 

1 In [HI [HI [20] the normalization of the scalar operator which breaks conformal invariance was such that 
its two point function vanishes in the limit A — > d/2. This also results in a vanishing speed of sound in this 
limit. The conventions in this work follow that of [2T]. To go from the conventions of [mESHD] ^° ^ ne ones 
in [21] one should divide the source and expectation value of the scalar operator by A — d/2. See section 2.2 
of [H] for details. 

2 In particular, one should make the identification A 2 = 2m 2 or A 4 = \m\ for the fermionic and bosonic 
deformations discussed in [2"2"] . 
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where P is the pressure, related to the energy density through the equation of state, r\ is the 
shear viscosity and ( is the bulk viscosity. The projection operator P^ u is given by 

pv» = jf + u »u v (8) 

and a^ u is given by 

= P m P vfi (d aU(3 + df,u a ) - -^P^P^d a up . (9) 

Greek indices run from to d — 1 and are raised and lowered with the Minkowski metric 
r\ = diag(— 1, 1, . . . , 1). By 0(d 2 ) I mean terms which contain two derivatives of the hydrody- 
namic variables, for example (J fJLa a a u . See [23] for a detailed explanation of higher derivative 
terms and the derivation of (17)). 

In a conformally invariant theory the stress-energy tensor is traceless, = , which 
implies e = (d — 1)P and ( = 0. In this note, conformal invariance will be broken by 
sourcing a relevant operator A of conformal dimension A, i.e, by adding to the boundary 
theory action a term — ^-j J A d ~ A 0&d d xQ Once the theory is not conformal, the trace of the 

2 

stress-energy tensor is given by 

TZ = -^4(0 A )A d - A + A. (10) 

The relation (flOj) can be derived by holographic means, as in |24| , or by computing the beta- 
function for the source term for Oa, as in [25]. For certain A there will be an anomalous 
contribution, A, to the trace of the stress tensor. Since A is independent of O it will not 
contribute to the speed of sound or bulk viscosity and will be ignored in what follows. 
Comparing (Tl0|) with ([7j) one finds that 

_ e + {d _ i )P _ {d _ i )CdaU <* + o{d 2 ) = -^^(0 A )A d ~ A , (11) 

^ 2 

so the bulk viscosity can be computed by extracting the term proportional to the gradient 
of the velocity field from the right hand side of ffTT]) . Such a computation will be carried 
out in the remainder of this work for a class of strongly coupled gauge-theory plasmas in the 
planar limit, which can be described using the gauge-gravity duality. 

3 As in [21], the factor (A — d/2)^ 1 has been introduced in order to have a non vanishing vacuum expec- 
tation value for OaOa in the A — * d/2 limit. 



3 



In detail, consider theories whose dual bulk action is given by 

S = i / ^ { R ~ {d ~ 1)d ~ + V{<f)) ) dd+±X (12) 

and where conformal invariance is softly broken by sourcing the operator Oa dual to the 
field 0. Working in a coordinate system where the line element near the asymptotically AdS 
boundary (located at r — ► oo) takes the form 

d-l 




-dv 2 + ^(dx*) 2 + 2dvdr , (13) 



the near boundary series expansion of the bulk field <ft is given by [21] 

- { ° A) -r"* (1 + O(r^)) + ^r*- d (l + O(r^)) +... (14) 



(2A-d) v ; A 
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where m 2 = 2V"(0) = A(d — A). The expectation value of the operator Oa dual to is 
given by (Oa), while A d ~ A is a source term for Oa- Once a solution to the equations of 
motion following from (1121) is found, the expectation value of Oa can be read off of (j!4p . 

A thermal state of the boundary theory corresponds to a black hole geometry [2^|. In 
the large temperature limit, one can keep only the leading terms in a A/T expansion of 
the solution to the equations of motion. The near boundary behavior of the scalar field <j) 
depends linearly on A d ~ A so to leading order one is allowed to keep only terms linear in 0. 
Thus, the Einstein equations are simply 

Rmn dg mn (15) 

where m, n — 0, . . . , d, and the equation of motion for the scalar is 

V 2 0-m 2 = O. (16) 

So the problem reduces to finding a black hole solution to the Einstein equations in the 
presence of a negative cosmological constant and finding a scalar field propagating in such 
a black hole background. It should be noted that this analysis can be easily generalized 
to several scalar fields: if the kinetic terms for the scalars are canonical then keeping only 
quadratic terms in the scalar potential amounts to decoupling the scalars from one another. 
To obtain the high temperature limit of the bulk viscosity of these theories I use (fTTl) 
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and compute the expectation value of Oa i n a black hole background which is dual to a 
configuration with d a u a 7^ 0. But before tackling this problem, it will be instructive to 
compute a simpler quantity, the speed of sound, c s . In a stationary configuration, where the 
spatial components of the velocity field vanish, one can use c 2 = dP/ de, to obtain 

<t) = - 71 ^ A^(Oa) (17) 



d-l 7 (d_i)(A-f) 

from ( fTTj) . In order to evaluate the right hand side of ( TT71) . one needs to compute the 
expectation value of the operator (Oa) dual to in a background dual to a stationary 
thermal state. Such a background is given by the AdS-Schwarzschild black hole solution 
whose line element is 

T 2 ( d-i \ 

ds 2 = — -f(r)dv 2 + ^{dx 1 ) 2 + 2dvdr (18) 

where 

and b = dj (4irT) with T the Hawking temperature of the black hole, and also the temperature 
of the thermal state in the boundary theory. The equation of motion for the scalar in this 
background is 



d r (r(l - (rb) d )d r (p) 



d-lJjd-2 



+ A(A-cf)0 = O. (20) 



The solutions to (120|) are 

<P = P_ 1+ a (-1 + 2(rb) d ) and <p = Q_ 1+ a (-1 + 2(rb) d ) (21) 



where P and Q are Legendre functions of the first and second kind. The solution to (1201) 
which is finite everywhere except at the black hole singularity and takes the asymptotic form 
(TTJJ near the boundary is 

<P(rb) = (A&)<*- A -^_ P _ i+# (_! + 2{rh y) . (22 ) 
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Expanding (122|) near the asymptotically AdS boundary and using (fT4"j) one can evaluate 

(0 A ) = [d) 2 l ^ (Ab) d ~ A . (23) 

Since the background is that of a conformal theory then to leading order in A/T the energy 
density e should scale as e ~ T d+1 . Using b = d/(4iiT) and (jTTl) . one obtains ([3]). 

Now I turn to the slightly more difficult problem of computing the bulk viscosity. As 
stated earlier, instead of the traditional approach to computing the bulk viscosity via the 
Kubo formula as was done in, for example, [23 [221 [28l El [291 [20], I will use (TTT1) in a 
background with non vanishing d a u a . Such black hole backgrounds were constructed in 
[301 I3T1 132"! [33] : omitting corrections of order 0(d 2 ), they are given by the line element 

ds 2 = ds 2 Q) + ds 2 1} (24a) 



where 



with 



and 



ds 2 0) = -2u^dx^dr - r 2 f(br)u^u v dx^dx u + r 2 P^ v dx^dx v (24b) 



ds 2 {l) = (^2r 2 (r^bF(br) + ^rd a u a u^u u - ru a d a v? {P^u„ + P lu u^j dx»dx v (24d) 



with 



F(br) = - m(fcr) + V ln ( fer X -K \ (24e) 
In f )24ep . the x n are the d — 1 roots of the polynomial 5^^=o^ n , and 6 and must satisfy 

a, (&-' {dux + 5/) - r^V) = o{d 2 ). (25) 

Note that b and u M now depend on the transverse coordinates x M . Put differently, (12~4"|) is a 
solution to (fl5l) to order 0(d 2 ) (as long as the boundary theory stress tensor is conserved, 
as in ([25]) ) and describes a configuration of the boundary theory with a non trivial flow. See 
[34] for a review. 

With the 0(9) solution to ([T5]) at hand one is left with solving (TIB]) in the background 
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(|24|) . Similar to the notation in f)24ap this solution can be written in the form 



+ , (26) 

where (p^ solves the scalar equation up to 0(d). Since at this order the metric is nothing 
but a boosted version of (JTBl . (j)^ is given by equation (1221) . 



2T 

4>V>(rb) = (A&) d " A ^^yP_ 1+ A (-1 + 2{rb) d ) . (27) 



Inserting ( |26|) and ( 1271) into ( |T6i) I find that 0^ must satisfy 
d r (r(l - (r6) d )9 r 0W) 



\d) r -(d-i)/2 dr ( r (i-i)/2p_ i+A (_i + 2(r6) d )) u a <9 Q ((A6) d " A ) (28) 

where ( 125]) was used to trade the gradient of the velocity field with the directional derivative 
of the energy density. After some guesswork, one finds that 

(fp-\rb) = R(rb)P^ 1+ A (-1 + 2(rb) d ) + CQ_ 1+ a (-1 + 2{rb) d+1 ) (29) 

with 

2r(-l 2 r°° v d ~ 2 

R(x) = -^j^-bu«d a ((Ab) d - A ) J x ^r^dy (30) 

and C a constant is a solution to (128]) which vanishes at the asymptotically AdS boundary. 
The function R(x) can be written explicitly in terms of logarithms. Requiring to be 
finite at r = 1/b implies 

C = -^P^L j bu«d a {(Ab) d ~ A ). (31) 

Expanding (I26I) as in (fT4l) one can extract the gradient corrections to (0&). Inserting this 
expression into (ITTj) and using (12"5|) one obtains (jl]). 

One would have hoped that using this method the second order transport coefficients 
of the nonconformal theory (see, for example, [35]) could be computed analytically. Unfor- 
tunately, it is difficult to solve the equations of motion for the 0(d 2 ) corrections to <fi. A 
similar problem arises if one tries to compute the bulk viscosity at finite chemical potential 
from the solutions in [SHI E3 EE] • 
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